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Abstract 



Using finite element software developed for metal cutting by Third Wave Systems we 
investigate the forces involved in chatter, a self-sustained oscillation of the cutting tool. 
The phenomena is decomposed into a vibrating tool cutting a flat surface work piece, and 
motionless tool cutting a work piece with a wavy surface. While cutting the wavy surface, 
the shearplane was seen to oscillate in advance of the oscillation of the depth of cut, as were 
the cutting, thrust, and shear plane forces. The vibrating tool was used to investigate process 
damping through the interaction of the relief face of the tool and the workpiece. Crushing 
forces are isolated and compared to the contact length between the tool and workpiece. We 
found that the wavelength dependence of the forces depended on the relative size of the 
wavelength to the length of the relief face of the tool. The results indicate that the damping 
force from crushing will be proportional to the cutting speed for short tools, and inversely 
proportional for long tools. 

keywords: metal cutting, chatter, process damping, shear plane dynamics. 

AMS classification numbers: 74S05,74H45,74H15. 
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1 Introduction 



In this paper wc investigate the forces produced in dynamic metal cutting, where either the 
chip load varies in time, or the tool itself oscillates. Chatter is a consequence of an instabil- 
ity in cutting that leads to self-sustained vibrations of the tool, hence it is composed of both 
an oscillating tool and an undulating workplace surface. We separate chatter into these two 
pieces by using finite element software for metal cutting developed by Third Wave Systems, 
Inc. called AdvantEdge, to simulate the processes. AdvantEdge is a validated software pack- 
age that integrates advanced dynamic, thermo-mechanically coupled finite element numerics 
and material modeling appropriate for machining processes. The simulation software provides 
accurate estimates of thermo-mechanical properties of the machining process such as cutting 
forces, chip morphology, machined surface residual stresses, and temperature behavior of the 
tool and workpiece. 

Recently Third Wave Systems has developed two new modules, one with custom work piece 
geometry and another with a vibrating tool, that increase the range of cutting conditions that 
can be simulated. The custom work piece geometry module allows for the design of arbitrary 

cutting surfaces in addition to the traditional flat surface. We study the effect of a varying chip 
load through a cut with this module, observing interesting shear plane dynamics and varying 
thrust and cutting forces. With the vibrating tool module, the cutting tool can be horizontally 
or vertically vibrated at a user specified frequency and amplitude. This allows us to investigate 
the effect of contact between the material and the relief face of the cutter, which generates a 
kind of process damping. 

In section 2 we examine the forces involved in cutting a workpiece with a sinusoidally 
varying surface and, as a consequense, a sinusoidally varying depth of cut. In section 3 we use 
the vibrating tool module to investigate process damping on the relief face of the tool during 
chatter. In section 4 we summarize our findings and suggest avenues for further research. 

2 Wavy Surface Runs 

The work piece geometry we use has a sinusoidal top surface which was generated by replacing 
the horizontal flat surface of the standard geometry with a sine function. The wavelength of 
the sinusoid was based on scaled parameters from experimental chatter measurements [8]. To 
ensure the feed was always greater than zero, the feed was set to be larger than the peak to 
valley amplitude of the sinusoid (3 mil). All these AdvantEdge runs were performed on AL 
7050 with a carbide tool with a 10 degree rake angle, and cutting edge radius of 0.7874 mil and 
cutting speed of 2680 SFM. Altogether three sets of runs were performed with three different 
feeds: 6,9, and 12 mil. 

From the animations of the Third Wave runs we saw that the shear plane angle during the 
cut oscillates as the tool moves through sinusoidally varying depth of cut (see figure 1). The 
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shear plane, denoted by the band of higher heat rates extending from the tool tip, is seen to lift 
up smoothly and shorten as the tool anticipates the thickest part of the chip. After reaching 
its maximum slope, the band snaps down and elongates to connect the tool tip with the point 
ahead on the surface that has the smallest chip thickness. The shear plane continues to extend 
between the tool tip and this minimum until that point, now a part of the chip, is crushed into 
the point on the surface with maximum chip thickness. The cycle repeats as the tool continues 
to encounter the oscillating chip thickness. The experiments we performed were designed to 
explore the effect of variation in depth of cut, shear plane angle, and shear plane length on the 
instantaneous cutting forces. 

2.1 Wavy Surface Force Predictions 

Our first observation is that the measured forces in x (cutting) and y (thrust), while oscillating 
at the same frequency as the chip thickness, lead it in phase. To illustrate this we plot the 
thrust and cutting force through the cut, overlaid with a sine wave of the same frequency and 
phase as the chip thickness, see figures 2 and 3. (The amplitude and DC offset of the sine wave 
were scaled to match the force for illustration). 

The phase lead makes sense after analyzing individual frames in the simulations. (See figure 
1.) The chip begins to bunch up and hence offer more resistance, which creates larger cutting 
forces prior to the maximum in the chip thickness. The forces drop before the tool tip encounters 
the next minimum of the chip thickness, in reaction to the shearplane snapping down to meet 
the minimum, which causes the chip to narrow down and offer less resistance. 

This phase lead of the force over the chip thickness has been observed in experiments carried 
out in Dr. Philip Bayly's lab at Washington University (for representative publications from 
Dr. Bayly's group see [2, 3, 4]). A tube cutting experiment with a vibrating tool was set up on 
a lathe, creating an essentially single degree of freedom system with one dominant mode. Tests 
were performed to determine the frequency response (FRF) of the cutting system, and these 
FRF were used to parameterize a single degree of freedom model for the cutting process. This 
model was developed in part by Brian Whitehead, who describes in his thesis [1] an attempt to 
include friction on the rake face of the tool by incorporating a weighting function that models 
the variation of forces as the chip slides up the rake face of the tool, after Stepan [5] . One of the 
parameters of the model is the contact length, L, between the tool relief face and chip. With 
positive L the weighting functions, once incorporated into the frequency domain model, act as 
low-pass filters, and induce a lag in the force behind the uncut chip thickness. Upon fitting 
this model with the tube cutting experiment FRF, they determined that the best fit came with 
using a negative L value, essentially making the forces on the tool anticipate the chip thickness. 
This does not make sense if L is considered to be the contact length on the rake face, but with 
L negative the force is being integrated over a region out in front of the tool, which induces the 
phase lead of the force over the uncut chip thickness. We believe that this is a manifestation 
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Figure 1: Sequence of snapshots from a Third Wave run, time increases to the right and down. 
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Thrust force for 12 mil cut, overlaid with chip variation. 
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Figure 3: Thrust force vs. distance into cut, overlaid with chip thickness (dashed hne). 
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Figure 4: Diagram of angles and forces for orthogonal cutting. 

of the shearplane oscillation effect we see in the Third Wave cutting simulations. 

It is plausible that while the forces lead the chip thickness, they may be in phase with the 
shearplane length, following a model popularized by Ernst and Merchant [6] , that assumes that 
the forces will be proportional to the shear plane length. The Merchant model predicts that 
specifically the shear plane force (Fg) is linearly dependent on the shear plane area, so Fg = twI, 
where r is the shear stress involved in cutting the material, w is the cutting width, and / is the 
length of the shear plane, which is oriented at the angle (p (refer to figure 4). Also from the 
figure note that Fg = ||R|| cos /3, where R is the total force on the tool. To determine /? wc rely 
on the traditional decomposition of R into vertical and horizontal components, Fy and Fx, the 
radial or thrust, and the tangential or cutting, forces respectively. If the angle that R makes 
with the horizontal is given by V = arctan(-pf^), then /3 = cj) + ip. The AdvantEdge program 
computes F^ and Fy, which also gives us a measure of the total cutting force as the tool moves 
through the cut, since ||R|| = ^F^ + F^ , as well as the angle tp. We can measure the shear 
plane length I, and the angle 0, for each step in the run, and since Fg = ||R|| cos(^ + we 
can test the Merchant prediction: 

Fs = \\K\\cos{(t) + tlj) =Twl. (1) 

The shearplane length was measured directly from the animation frames, where visualiza- 
tions of the heat production are best for identifying the shearplane region. This is subject to 
obvious error through the identification of a line that could be called the shearplane. We chose 
to measure a line that started at the tip of the tool and ended at the surface of the workpiece, 
parallel with the contours of heat production. See figure 5 for an example. The shear plane 
length oscillates along with the undulations in the surface of the workpiece, though not neces- 
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X(mm) 

Figure 5: How the shearplane angle and length is computed: the angle of the line shown here 
with the horizontal and its length are measured directly from the frames of the Third Wave 
animation. 

sarily with the same phase. Following this observation, the shearplane length as the tool moves 
through the workpiece was least squares fit to a sine wave with the same wavelength as the 
wavy surface. The result of this exercise for the 12 mil feed run is plotted in figure 6, and the 
fitted expression is 

l{x) = 4.6 + l.l(.151cos(27rx/40) +0.647 sin(27rx/40)). 

The units on shearplane length are arbitrary, and vary with the three runs examined in this 
section. The important feature is the phase of the oscillation. 

The shearplane angle can be measured and fit in a similar way, the result is presented in 
figure 7, and the fitted expression is 

(l)(x) = 0.49 + 0.1 cos (27rx/40) - 0.09sin(27rx/40). 

Similar fits were performed for the 6 and 9 mil cuts, and the units are radians. We also observe 
that the shearplane length oscillation leads the chip thickness oscillation by a noticeable amount, 
see figure 8. 

We compute ||R|| = y^F^ + Fy and /? = (/) + arctan(^) for each of the three runs as a 
function of distance into the cut (x). Expression (1) implies that l{x) ~ ||R(x)|| cos/3(x), so 
we conducted a least squares fit to ||R|| cos /3(x) with the expression ci + C2l{x), the result is 
presented in figure 9. The shear plane force is seen to lag the shear plane length, which fits 
with the observation that while the shearplane force leads the chip thickness, the shearplane 
length leads the chip thickness by an even greater amount, leading to the lag in shearplane force 
with respect to the shearplane length. Also note that we fit the force Fg not with cil{x), as is 
specified by Ernst and Merchant, instead it is a linear fit offset by a constant. This indicates 
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Shear Plane Length, 12 mil cut 




distance into cut (mils) 

Figure 6: Shear plane length measured from experimental snapshots, overlaid with fitted sine 
wave. 



Shear Plane Angle, 12 mil cut 
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Figure 7: Shear plane angle measured from experimental snapshots, overlaid with fitted sine 
wave. 
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Shear plane length compared to chip thickness, 12 mil cut 
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Figure 8: l{x) compared to chip thickness through the cut for 12 mil cut. The magnitude 
of the chip thickness is in mils, the shearplane length has been scaled and shifted vertically 
for comparison. The critical feature is the phase shift: the shearplane length leads the chip 
thickness. 

that there must be a nonlinear region near / = 0, a roll-off of the function so that Fs = when 
I = 0. This by itself does not contradict the Merchant model, which was meant to apply only 
in steady cutting finite chip thickness regimes. However, the Merchant formula does fail to 
completely capture the behavior of the forces in this dynamic situation. 

This procedure suggests plotting the data ||R(x)|| cos parametrically against l{x) to 
determine if this linear fit is appropriate, see figure 10. In this format it is clear that the 
shearplane force depends on something in addition to /, since it is a multi- valued. The fitted 
line is drawn for comparison for each of the three runs. 

The direct measured forces, and Fy can also be compared with the shearplane length 
oscillation during the run. The result is presented in figure 11, for the 12 mil chip thickness run. 
The cutting force also lags the shearplane length, similar to Fg, predictably since the cutting 
force is the largest component of the total force. The thrust force, however, comes very close 
to matching the oscillation of the shearplane length, a fact that could be explored further. 

To conclude this section, we would like to comment on an idea presented in Wellbourne and 
Smith [7] . They computed an expression for the cutting force involved in cutting a wavy surface 
by assuming that the force was proportional to the shear plane length also. However, they were 
limited to estimating the shearplane length by computing the length of a line, oriented at a 
constant slope, as it intercepts a sine wave (see figure 12). The length can be found by solving 
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Figure 9: Fs compared with fit to shearplane length for 12 mil, 9 mil and 6 mil depth of cut 
(D.O.C.) runs, the length is the smooth curve. 



11 



6 mil feed 



9 mil feed 




2 4 6 8 
shear plane length 

12 mil feed 



10 





1.5 2 2.5 

shear plane length 



3 4 5 6 

shear plane length 

Figure 10: ||R(x)|| cos/5(x) plotted parameterically with x against I, 6,9 and 12 mil feed. The 
shear plane length is in different units in each figure, the measured shearplane force is in pounds. 
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Cutting and thrust forces compared with shearplane length, 12 mil D.O.C. 
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Figure 11: and Fy plotted with I (the smooth curve), 12 mil feed. 
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Figure 12: How the shear plane length might change through a wavy surface cut. 
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shearplane length through cut, varying shearplane angle 
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Figure 13: Plotting the sliear plane length for varying shear plane angle ((p), through a wavy sur- 
face cut. The function plotted with a dashed line is the workpiece surface itself, for comparison 
with the chip thickness. 



a transcendental equation (the intersection of the moving line and a sine wave) numerically 
for varying shear plane angles (figure 13). From this you can see that the length leads the 
oscillation of the chip load, how much depending on the shear plane angle. The oscillation also 
becomes more triangular as the angle is decreased, it will have a discontinuity when the multiple 
roots to the transcendental equation appear, that is, for small enough shearplane angles. 

These effects certainly play into the form of the forces determined with experimentally 
measured shear plane lengths. But there are several confounding factors: the shear plane angle 
is not a constant, and more importantly, the surface itself is deformed during a cut, and the 
simple intersection of a line and a sine wave does not give the correct shear plane length, as we 
saw in the animations, see figure 1. 

3 Analyzing Crushing Forces 

As described in the Introduction, the new Third Wave module that allows the tool to oscillate 
either horizontally or vertically during a cut gave us an opportunity to explore the forces 
involved when the relief face of the cutting tool impacts the workpiece. In our simulations to 
investigate relief face crushing the cutting tool was vibrated sinusoidally in the vertical direction 
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Figure 14: Third Wave run (AL7050) with rehef angle of 6 degrees. 



with a 3 mil amplitude and varying frequencies. Two simulations were performed to calculate 
the crushing force for each set of cutting parameters, identical in all ways except for the relief 
angle of the tool. In one simulation, the relief angle was small enough that the entire relief face 
came into contact with the workpiece, causing the crushing of the workpiccc at the relief face 
of the tool. In the other simulation, the relief angle was made large enough so that there was 
no contact between the relief face of the tool and the work piece. 

To illustrate, a snapshot of a crushing run is shown in figure 14, where the tool has come 
into contact with the workpiece during the downstroke of the tool and left a flattened region in 
its otherwise sinusoidal path. This should be compared with a run with the relief angle drawn 
up from 6 degrees to 25 degrees, figure 15, where the tool tip leaves a sinusoidal workpiece 
surface in its wake. Assuming the forces involved in material removal and crushing combine 
linearly, the crushing force can be resolved by subtracting the forces from two such runs. The 
exercise is illustrated in figure 16, where the material being cut is AL7050, with a tool with rake 
angle of 10 degrees and cutting edge radius of 0.5 mils. In the two runs the wavelength of the 
oscillation is 80 mils, corresponding to a cutting speed of 2680 SFM and vibration frequency of 
6.7 kHz. The nominal feed is 1 mil, and vibration amplitude is 3 mil. The spatial wavelength 
A of the undulating cut is equal to the ratio of the vibration frequency v to the cutting speed 
V, that is A = ^. We focused our investigations on the vertical or thrust force in what follows. 
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Third Wave AdvantEdge 
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Figure 15: Third Wave run (AL7050) with rehef angle of 25 degrees. 

Wc found that the rehef edge length, relative to the wavelength of the oscillation, was critical 
in determining the behavior of the crushing force in these runs. Therefore we ran simulations 
with tool relief lengths of 10 mil, 30 mil and 100 mil, which were much smaller, somewhat 
smaller, and much larger than the range of wavelengths tested. The results are summarized 
below. 



3.1 Crushing force simulations 

Once again the AdvantEdge runs were performed on AL 7050 with a carbide tool with a 10 
degree rake angle and cutting edge radius of 0.7874 mils. The cutting speed was chosen to be 
2680 SFM, and vibration frequencies near 13.4 KHz, and amplitudes of 3 mils, typical of an 
axial-torsional vibration mode for a twist drill ([8]). In total we did 24 runs: 12 crushing, 12 
no crushing. For each length of tool relief face (10, 30 and 100 mils) we ran four wavelengths: 
40, 60, 80 and 100 mil. We found plotting the crushing force vs. the vertical amplitude of 
the tool during a run a useful way of visualizing the data. In figures 17-19 these results are 
presented for each tool relief length, with the varying wavelength runs plotted together for 
comparison. The general form of each plot is the same, during the downstroke of the tool, 
as expected, the crushing force increases up the angled portion of the loop, i.e. the loops are 
traversed counterclockwise with increasing time. Near the end of the downstroke the force 
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Figure 16: Computing the crushing forces from two Third Wave runs. 
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Long tool relief length (100 mil) 
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Figure 17: Thrust crushing force vs. vertical position of tool, long tool relief length, multiple 
wavelengths. 

begins to decrease, and drops to zero abruptly as the tool begins its trip back up to the top 
of the oscillation. During this phase the force is zero and the portion of the loop along the 
displacement axis is traversed. Immediately after the top of the oscillation is reached and the 
downstroke begins, the force starts to increase again, and another cycle is initiated. 

The three sets show a qualitative difference in shape, and a quantitative difference in maxi- 
mum force achieved relative to wavelength of the oscillation. The long relief length tool demon- 
strates the most easily understood force loop: the force increases approximately linearly as the 
tool moves into the material, reaching a maximum close to the lowest point of the traverse. 
In the short relief length tool runs the force increases approximately linearly until about one 
third of the way down, when it flattens out. The short tool relief length is saturated at this 
point: is has come into contact with as much material as it can, it cannot impact more length 
as it moves down. The intermediate length tool runs have a shorter flat region that depends 
on wavelength, with longer wavelengths showing some flattening. 

The quantitative difference between the three sets of runs concerns the dependence of the 
maximum force value vs. wavelength of oscillation. In Whitehead [1], a model for process 
damping due to rubbing on the relief face is proposed that depends inversely on the wavelength 
of the vibration, based on a suggestion by Tobias [9]. This translates to a term in his linear 
model that is directly proportional vibration velocity (which depends directly on vibration 
frequency) and inversely proportional to cutting speed. In figure 20 we plot the maximum 
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Intermediate tool relief length (30 mil) 




Figure 18: Thrust crushing force vs. vertical position of tool, intermediate tool relief length, 
multiple wavelengths. 



Short tool relief length (10 mil) 
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Figure 19: Thrust crushing force vs. vertical position of tool, short tool relief length, multiple 
wavelengths. 
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wavelength 

Figure 20: Maximum crushing force vs. wavelength of oscillation. 

force for each run vs. wavelength in the three cases. Though the data is too sparse to fit 
any detailed model, it is clear that the force is inversely proportional to the wavelength only 
in the short relief length runs, the long relief length runs demonstrate a direct dependence on 
wavelength, and in the intermediate length runs the maximum force appears to be indifferent 
to wavelength. The dependence of the damping on cutting speed and vibration frequency thus 
will depend on the length of the relief face relative to the vibration wavelength, and given that 
chatter vibrations are mainly high frequency, the long tool data will be the most relevant for 
estimating chatter process damping. 

3.2 Crushing force modeling 

The magnitude of the crushing force will most certainly depend on the amount of contact 
between the relief face of the tool and the material. We can analytically model the variation of 
the contact during a crushing thrust of the tool by determining the contact length, as it depends 
on frequency of the oscillation, height of oscillation and the cutting speed. The contact length 
at each point along the cut can be calculated from the intersection of a tool edge with the 
work piece (see figure 21). The sinusoid in the figure is the path that the tool takes during an 
oscillation, the chip is not shown in the diagram. A small program in MATLAB calculated this 
intersection length in the figures that follow. 

The contact length during an oscillating cut is plotted vs. vertical position of the tool tip, 
see figure 22. This representation can be directly compared to the force vs. vertical position 
plot for the 3rd Wave run with a 100 mil relief length tool, figure 17. The contact length and the 
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Figure 21: Illustration of contact length during crushing, long tool relief length. 



force have similar dependence on wavelength, both increase with increasing wavelength at the 
bottom of the downstroke of the tool (the highest point on the loop). There is also a "roll-over" 
at the top of the downstroke, where the shorter wavelength force/contact length docs not begin 
to grow as soon as the longer wavelength case. This is a consequence of the relative slope of 
the tool path to the relief angle of the tool. A similar diagram for the short relief length (10 
mil) tool is shown in figure 24. The contact and the crushing forces in this case saturate when 
the tool relief length is fully engaged with the workpiece, but an obvious difference between 
the two is the invariance the maximum contact length shows with respect to wavelength. The 
maximum forces, however, decrease with increasing wavelength for the short tool. 

With this evidence in mind we can attempt to fit the crushing force data from the experi- 
ments with a model that has a linear dependence on contact length. If X is the contact length, 
Y = mX + 6 is the crushing force, and m and h are fitted via least squares. The results depend 
on the wavelength of the oscillation and the relief length of the tool, and are tabulated below. 

Long relief length tool (100 mil): 



wavelength (mils) 


fitted line 


error (lbs) 


40 


Y = 1.102 


* 


X + 0.397 


2.34 


60 


Y = 1.132 


* 


X - 1.539 


3.75 


80 


Y = 1.184 


* 


X - 9.698 


3.88 


100 


Y = 1.305 


* 


X - 20.493 


3.74 



Intermediate relief length tool (30 mil): 
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Contact Length vs. position of tool, long relief length 
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Figure 22: Contact length during tool oscillation, long tool relief length (100 mils). 



Contact length vs. position of tool, intermediate relief length (I 
35 1 1 1 1 1 — I 1 




vertical position of tool (mils) 



Figure 23: Contact length during tool oscillation, long tool relief length (30 mils). 
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Contact length vs. position of tool, short relief length 
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Figure 24: Contact length during tool oscillation, short tool relief length (10 mils). 



wavelength (mils) 


fitted line 


error (lbs) 


40 


Y = 1.002 


* 


X + 3.628 


3.02 


60 


Y = 1.431 


* 


X - 6.472 


2.54 


80 


Y = 0.963 


* 


X - 7.128 


2.47 


100 


Y = 0.977 


* 


X - 7.877 


2.62 



Short relief length tool (30 mil) : 



wavelength (mils) 


fitted line 


error (lbs) 


40 


Y = 0.548 * X + 1.501 


0.84 


60 


Y = 0.234 * X + 1.350 


0.51 


80 


Y = 0.207 * X + 1.484 


0.12 


100 


Y = 0.167 * X + 0.957 


0.33 



The slope constant in the linear relationship gives the scale factor between force (in lbs) 
and contact length (in mils) and thus has units lbs/mil of contact length. This conversion 
factor varies somewhat within runs (mean and standard deviation are 1.18 it 0.089, 1.093 it 
0.226,0.289 =b 0.175, for the long, intermediate and short relief length runs, respectively) The 
y intercept is determined by nonlinearity in the data at the top and bottom of the cycle in the 
data, which causes the linear region to move up or down. 

To illustrate these results we plot the scaled contact length vs. vertical position of the tool 
along with the crushing force for a comparable run, for the three different relief lengths and 
varying oscillation wavelengths (figures 22,23, and 24). The long tool force loop shows a close fit 
with contact length, except near the top of the oscillation when the contact length grows more 
quickly from zero. The intermediate tool shows a better fit with data at top of the oscillation, 
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Figure 25: Crushing force and contact length during tool oscillation, long tool relief length, 80 
mil wavelength. 

though the force does not flatten out in the same way as the contact length. The short tool 
forces are undercut the contact length, the contact length grows more quickly at the top of 
the oscillation, again, and the forces begin to decrease from the plateau before the contact 
length comes down. All these point to additional factors coming into play in the prediction of 
crushing forces. The assumption of linear dependence on contact length holds primarily for the 
mid-region of the downstroke of the tool, excluding the lift-off at the bottom and turn-around 
point at the top. 

To illustrate these results we plot the scaled contact length vs. vertical position of the 
tool along with the crushing force for a comparable run, for the three different relief lengths 
and varying oscillation wavelengths (figures 25, 26, and 27). The long tool force loop shows 
a close fit with contact length, except near the top of the oscillation when the contact length 
grows more quickly from zero. The intermediate tool shows a better fit with data at top of 
the oscillation, though the force does not flatten out in the same way as the contact length. 
The short tool forces are undercut the contact length, the contact length grows more quickly 
at the top of the oscillation, again, and the forces begin to decrease from the plateau before 
the contact length comes down. All these point to additional factors coming into play in the 
prediction of crushing forces. The assumption of linear dependence on contact length holds 



24 



Intermediate relief tool - 30 mil, 80 mil wavelength 

40 I , , , , , 




o 
O 

-5 I ' ' ' ' ' 

-3-2-10123 

Vertical Position of tool (mils) 



Figure 26: Crushing force and contact length during tool oscillation, intermediate tool relief 
length. 
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Figure 27: Crushing force and contact length during tool oscillation, short tool relief length. 
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primarily for the mid-region of the downstroke of the tool, excluding the lift-off at the bottom 
and turn-around point at the top. 

4 Discussion and Conclusions 

Our investigations with Third Wave simulations were designed to isolate different effects present 
during regenerative chatter vibrations in metal cutting: the vibration of the tool and the 
oscillatory variation in the chip thickness. The tool vibration feature of the package was used 
to investigate process damping due to crushing of the workpiece material on the relief face of 
the tool, and the custom workpiece feature allowed us to set up a sinusoidally varying chip- 
load. Simulations with the later demonstrated an oscillation of the shear plane as the tool 
progressed, and we found that the cutting, thrust and shear plane forces all were in advance 
of the chip thickness. The shear plane force, proposed to be directly dependent on shear plane 
length by Ernst and Merchant [6] , was found to lead the shear plane length as well, as the forces 
anticipated the decrease in chip load. This anticipation is due primarily to the deformation of 
the chip as it encounters the tool, not surprizingly the process is more complex than predicted 
during dynamic cutting. 

An important cause of process damping is the interaction of the relief face of the cutter 
with the workpiece material. Here we studied the thrust forces involved in such damping by 
running an oscillating tool in AdvantEdge for a variety of oscillation wavelengths and tool relief 
lengths. We found that the wavelength dependence of the force was determined by the length 
of the relief face of the cutter relative to the wavelength of the oscillation: for long relative relief 
length the maximum force increased with increasing wavelength, for short relative relief length 
it decreased, and for intermediate relative relief length the dependence was indeterminent. 
This expands the intuition of Whitehead et al., [1] who used a process damping model that 
depended inversely on wavelength, as it would for a short tool. Since the wavelength of the 
oscillation equals the cutting speed divided by the oscillation frequency, this means the damping 
is proportional to the cutting speed for short tools, and inversely proportional to the cutting 
speed for long tools. This, however, assumes that the frequency of oscillation is independent of 
cutting speed, which is most likely not the case, given the form of the linear stability diagrams 
for regenerative chatter (for an example related to this study see [8]). 

We also found that the force depended linearly on contact length between the relief face of 
the tool and the workpiece material during the downward stroke of the tool. Other factors come 
into play during the initial moments of the downward stroke, and through the "lift-off of the 
tool from the trough of the oscillation. This suggests that this form of process damping can be 
explained purely geometrically, where the relevant factors are oscillation amplitude and relative 
size of relief length to wavelength. Incorporating this into a model that would determine 
stability of the steady cutting state to oscillations would not be straight-forward, however. 
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since the damping itself depends on the wavelength of the instability induced (typically the 
wavelength is determined by the physical parameters of the system, not a priori) . 

In summary, this research highlights the importance of dynamic and nonlinear effects on 
instabilities in metal cutting, even in the case of orthogonal cutting. Forces that act in anticipa- 
tion of chipload increases and process damping that depends in a complicated way on vibration 
frequency are examples of effects that could be studied further, both in the laboratory and with 
large scale and predictive computer modeling. 
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